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Abstract 

Possible variations in the dynamical behaviour of LiYF4 due to its several struc- 
tural changes under pressure are examined by making use of the complementary 
techniques of quasi-harmonic lattice dynamics and molecular dynamics simulations. 
The phonon spectra in the entire Brillouin zone together with the respective Gibbs 
free energies are calculated for the three high-pressure polymorphs of L1YF4 (that 
are stable at T = 0) with a view to better understand their relative stabilities as 
functions of pressure and temperature. The present work predicts anomalous ther- 
mal expansion at low temperatures in phases I and Ha while irreversibilty of phase 
II — > phase III transition on subsequent pressure release. Molecular dynamics sim- 
Q\ ■ ulations provide qualitative impressions about a temperature-driven second-order 

<^ transition and also of kinetic effects in the subsequent pressure-driven first-order 

I phase transformation. 
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1 Introduction 

O . 

O ■ Ternary rare-earth halide LiYF 4 , which has drawn considerable interest[l,2,3,4,5,6,7,8,9,10] 

as a laser host material, displays several structural phase transitions[4,6,8] 
upon compression. While isostructural CaWC>4 has a scheelite (IAi/a) to wol- 
framite (P2i/c) kind of high-pressure phase transformation, L1YF4 displays 
an intermediate fergusonite phase (viz. 12/ a, above room temperature and 
P2i/c, in the low temperature region) associated with a soft phonon mode[5]. 
Further, the axial ratio in both the compounds evolves in a different way un- 
der the effect of pressure [10]. Recently, the post-fergusonite phase of LiYF 4 
has been predicted as wolframitelike monoclinic one[8,9]. 



> 



c3 



1 Present address: Department of Physics, Indian Institute of Technology, Kanpur 
208016, India; Electronic address: asen@iitk.ac.in 



Preprint submitted to Elsevier Science 



2 February 2008 



With the advent of high-energy synchrotron X-rays and sophisticated 
pressure cells, high-pressure lattice dynamics has gained considerable momen- 
tum in recent years as it enables a direct comparison of experimental findings 
obtained either by the inelastic scattering processes or through some other 
techniques like ultrafast X-ray diffraction [11]. In this work, we have carried 
out extensive lattice dynamical calculations in the quasiharmonic approxima- 
tion (QHA) for various individual high-pressure phases (e. g. I, Ha and III in 
Fig. 1) of LiYF 4 that may help analyze experimental observations. The same 
interatomic potential as developed earlier [3] for lattice dynamical calculations 
under ambient conditions is made use of here again for high-pressure studies. 
Gibbs free energies are calculated for various high-pressure phases of LiYF4 
in order to have a better understanding of their relative stabilities in terms 
of volume compression and vibrational entropy. We further examine their re- 
spective phonon dispersions and phonon density of states. The second-order 
phase transitions and the kinetic effects in the first-order phase transforma- 
tions are specifically investigated by molecular dynamics simulations. Finally, 
a pressure-temperature (P-T) phase diagram associated with LiYF4 is plot- 
ted to qualitatively gauge the underlying physics of phase transitions at the 
microscopic level. 



2 Vibrational Properties 

2.1 Raman and infrared active modes 

From the previous works[2,3,5], we know that the ambient scheelite phase 
(74i/a, Z=4) of LiYF4 has a total of 36 vibrational degrees of freedom per 
primitive cell. However, the two high-pressure phases (Ha and III in Fig. 1) 
possess 72 phonon modes at each wavevector due to doubling of the primitive 
cell. Since both these high-pressure phases of LiYF 4 belong to the same space 
group(P2i/c, Z=4), the common irreducible representation of the phonons at 
the zone center is given by 

Y: 18A 9 +18A U +185 9 +185 U , 

of which A g and B g modes are Raman active while modes of ungerade sym- 
metry (viz. A u and B u ) remain infrared active. Although our calculated high- 
pressure phases have 36 Raman active modes in all, the ambient fergusonite 
structure has only 18 such frequencies. It may be further noted that the de- 
generacy of the Eg and E n modes which were present in phase I (scheelite) of 
LiYF4 is lifted in the high-pressure phases. Table 1 yields a comparison of our 
calculated zone centre normal modes between phases Ha and III at respec- 
tively 13 and 18 GPa. Phase Ha has a soft A g mode at 18 cm -1 which hardens 
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considerably to 92 cm^ 1 in phase III. Most other modes harden slightly in 
going from phase Ha to phase III. The vibrational frequencies for each of the 
phases are computed by diagonalization of the respective dynamical matrices 
using the software DISPR[13], developed at Trombay. 

2.2 Phonon dispersion relations 

To compare the phonon dispersion in various phases, we resort to the 
common high symmetry direction of A, which is labeled as the c axis for phase 
I (I4i/a) while it is the b axis for the other high-pressure(monoclinic) phases 
of similar space group symmetry (i. e. P2i/c). From our group theoretical 
calculations, we obtain for the irreducible representations of all the normal 
modes in the three phases as follows: 

Phase I -> A: 8A1 + 8A 2 + 10A 3 (A 3 being doubly degenerate) 
Phase Ha and Phase III -> A: 36Ai + 36A 2 

Phonon dispersions of various polymorphs in LiYF 4 are displayed in Figure 2 
along with the available inelastic neutron scattering data. [2] The experimental 
data [2] for the phase I at zero pressure fit well with the lattice dynamical 
calculations [3]. Due to doubling of the unit cell in going from phase I to phase 
Ha, the Brillouin zone in phase Ha is halved and so there is an apparent folding 
back of the dispersion branches from zone boundary (of phase I) to zone center 
(of phase Ha). This leads to several optic phonons at the zone center in phase 
Ha at low energies. We find a very low energy zone-center optic mode in 
phase Ha, indicating the onset of a possible high-pressure phase transition at 
a pressure close to 13 GPa. This espouses our previous MD observations [8]. 
The soft-alike phonon branch of phase Ha may be identified as the longitudinal 
A g mode [Table 1] as noted earlier. A crossover between an acoustic phonon 
mode and a low frequency optic mode is observed in phase Ha and also in phase 
III. We further notice that two acoustic phonon branches of A 2 symmetry in 
the phase III are too closely dispersed to remain distinguishable. 

2.3 Phonon density of states 

To have an overview of the range and extent of various phonon modes in all 
the three phases, we integrate over all the phonons with an energy resolution 
of 0.5 meV at each wave vector on a 16 x 16 x 16 mesh within the irreducible 
Brillouin zone. As Figure 3 suggests, the energy distribution extends by nearly 
10 % in switching from phase I to phase Ha while the change in the extent is 
negligible for the phase III transition. In order to check whether it is due to the 
pressure effect, we compare the respective frequency distribution for phase Ha 



3 



and phase III at the same pressure of 13 GPa. One may notice that though the 
distribution patterns slightly differ, the extent remain quite the same. If one 
is further interested in looking for individual contributions of the constituent 
atoms to the entire frequency distribution, partial density of states (PDOS) 
are to be estimated by atomic projections of the one-phonon eigenvectors. 

The labeling of the crystallographic axes viz. a, b, and c varies in high- 
pressure phases as per convention adopted for specific space groups. In view 
of these differences, we choose right-handed orthogonal x, y and z axes as 
follows: z along c and x along a in phase I (IAi/a); z along b and x along c in 
phases Ha and III (P2 1 /c). Our calculated PDOS for various atoms in various 
phases are portrayed in Figure 4. Y atoms are found to contribute largely up 
to 60 meV while Li atoms make it (especially with polarization along x and 
y) in the higher energy side of the density of states (DOS). However, F atom 
contributions are spread over the entire energy range as if to replicate the 
total DOS to a considerable degree. Such variations in atomic contributions 
are partly due to the mass effect. 



3 Thermodynamic Properties 

3. 1 Heat capacity and Debye Temperature 

From the DOS calculations derived out of our simulated data for all the 
three phases of LiYF 4 , we compare the respective heat capacities at constant 
pressure(Cp) as a function of temperature. Figure 5(a) demonstrates how Cp 
changes for different high-pressure phases. We observe that although initially 
constant pressure heat capacity for phase III runs lower, it goes up above 
the room temperature and in the process, even surpasses phase Ha. The rea- 
son is that as we go on increasing temperature (T), the difference Cp(T) - 
Cy(T) [=ay(T)BVT] (where Cy(T) is the constant volume heat capacity 
and B is the isothermal bulk modulus) becomes significantly large, as inset of 
Figure 5(a) suggests, due to the existence of much higher thermal expansion 
ay(T) in phase III than that in phase Ha. A detailed discussion on ay(T) 
follows in the next subsection. 

We compare the Debye temperature D (T) for the various phases [Fig- 
ure 5(b)] as derived from the calculated Cy(T). It is observed that in the very 
low temperature region (say, < 20 K), 9 a differs a lot from phase to phase 
(about 14 % increase for I— >IIa and about 26 % increase for Ha— TII) while 
at higher temperatures (say, 300 K), it comes closer for phases Ila and III as 
difference gets minimal (about 2 % only). 
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3.2 Griineisen parameter and thermal expansion 



Lattice anharmonicity, which leads to a volume dependence of phonon 
frequencies^), is described by the mode Griineisen parameter[14] given as 



It is the only kind of anharmonicity that can be taken care of within the frame- 
work of quasiharmonic approximation. Moreover, Griineisen parameter is an 
important quantity as it describes the thermoelastic behaviour of materials at 
high pressures and temperatures. It has both the microscopic and macroscopic 
definition, the former relating to the vibrational frequencies of atoms in a ma- 
terial (Eq. 1) while the latter, to familiar thermodynamic properties such as 
heat capacity and thermal expansion. Unfortunately, the experimental deter- 
mination of 7i, defined in either way, is often not so easy, since the microscopic 
definition requires a detailed knowledge of the phonon dispersion spectrum of a 
material, whereas the macroscopic definition calls for experimental measure- 
ments of thermodynamic properties at high pressures and temperatures. In 
this perspective, theoretical model calculation may be of great relief. Figure 6 
displays our calculated average 7, for various phonon energies and in various 
phases. It can be seen that below 10 meV, there is a significant number of 
modes with negative Griineisen parameter in phase I as well as phase Ila, 
while phase III modes have all positive 7*. 

The effect of pressure on the volume coefficient of thermal expansion(o;y) 
can be studied through mode Griineisen parameters (7^) in the entire Brillouin 
zone. In the quasiharmonic approximation, each phonon mode of energy Ej 
contributes to the thermal expansion by way of (^)7iCvj, where Cyi denotes 
the specific heat (at constant volume) contribution of the i th mode and V 
the lattice volume. This procedure suits well to the ambient ffuoroscheelite 
system because explicit anharmonicity (which arises mainly out of thermal 
amplitudes) is not very significant [5] compared to the implicit effect that in- 
volves phonon frequency change with volume (as obtained from 7,). 

Figure 7(a) demonstrates the variation of ay as a function of tempera- 
ture. It is also apparent that the scheelite as well as the initial high-pressure 
phase of LiYF 4 have negative thermal expansion in the low temperature limit 
(below 100 K). It may be because 7, has large negative values for phonons 
below 10 meV in these two phases. However, the third phase has no such 
anomaly, since it has all positive gammas for phonons of all energies (Fig. 6). 
The contributions to ay from phonons of different energies (corresponding to 
various phases of LiYF 4 ) are displayed in Figs. 7(b), (c) and (d) respectively. 



din iOi 




li = 



din V 



5 



We observe that at low temperatures (e.g. 20 K), contribution of modes upto 
10 meV are quite significant to ay, but as temperature is increased (e. g. 300 
K) higher energy modes get more populated and hence, contribute in a large 
way to the volume thermal expansion. 



3.3 Mean squared displacements and thermal anisotropy 



In order to gain some insight over how the phonons of various energies 
are polarized in various phases of L1YF 4 , we plot (see Figure 8) the partial 
contributions of these phonons to the mean squared thermal amplitude of 
the constituent atoms. The mean squared displacement of atom k along a 
direction is given by 

U aa (k,T) = (u 2 ka ) T = A— ] {n + 1/2) g ka (")dou, (2) 



exp (j?f) — 1 ; gka(y>) is the partial density of states associated 



where 

with an atom k whose mass is m^; A is the normalization constant such that 
/ gka{u)duo = 1. As Figure 8 suggests, the modes at very low energies involve 
equal displacements of all the atoms that correspond to the acoustic modes. 
Interestingly, between 2 and 9 meV, Y and F atoms have larger amplitudes 
than what relatively light-weight Li atoms possess. It may be noted that the 
basic structure (scheelite) of LiYF 4 comprises a pair of strongly bonded LiF 4 
tetrahedra and loosely bonded YF 8 polyhedra[10]. Larger amplitudes of F 
atoms in the first two phases mark the presence of librations of the LiF 4 
tetrahedra. 



Our calculated values of equivalent isotropic thermal parameters for Li, 
Y and F atoms in the ambient (i. e. P=0) scheelite phase of LiYF 4 are found 
to be 17, 9 and 14 x 10 _2 A 2 as against the respective experimental [15] values 
of 20, 10 and 17 x 10- 2 A 2 . A detailed comparison of the anisotropic ther- 
mal parameters (U aa ) at different temperatures among the different phases of 
LiYF 4 is further given in Table 2. It may be interesting to note that while Y 
and F atoms have comparable atomic displacements along the three directons 
in all the high-pressure phases of LiYF 4 , Li atoms, in contrast, show larger 
anisotropy along the z direction in phase I, along the y direction in phase II 
and again along the z direction in phase III. 
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3.4 Free energy and phase stability 



A pressure-temperature phase diagram generally portrays various equi- 
librium phases at constant temperature (T) and pressure (P) with the lowest 
Gibbs free energy (G). From phonon calculations, the temperature dependent 
vibrational free energy at various hydrostatic pressures for various phases of 
LiYF4 can be estimated in the quasiharmonic approximation. Thermodynam- 
ically, we may write[16] 

G n (P,T) = U n -TS n + PV n (3) 

where U n , S n and V n refer respectively to the internal energy, the vibrational 
entropy and the lattice volume of the n th phase. The lattice excitations are 
treated in this work within the quasiharmonic approximation where the full 
Hamiltonian at a given volume is approximated by a harmonic expansion 
about the equilibrium atomic positions and anharmonic effects are included 
through the implicit volume dependence of the vibrational frequencies. Hence, 
the vibrational Gibbs free energy (G" b ) is found to be satisfyingly accurate. 

To include vibrational effects in the present phase diagram, we have cal- 
culated dynamical matrices separately for each of the three phases (viz. I, Ha 
and III) at pressure intervals of 2 GPa on a 4 x 4 x 4 mesh throughout the 
irreducible Brillouin zone comprising 64 wavevectors. In order to rationalize 
the behaviour of this dynamical simulation, we calculate [17] the enthalpy vs. 
pressure curves for all the three known structures of LiYF4 and notice that 
the enthalpy changes (due to internal energy and volume) are predominant in 
pressure- driven transitions over free energy changes (due to vibrational energy 
and entropy). Perhaps this is why most of the first-principles phase diagram 
calculations of pressure-driven phase transitions [18] do not include vibrational 
entropies, though configurational contributions are always taken into account. 
The other reason for this may be that vibrational entropy differences between 
phases are assumed to be quite small as we come across in this particular case 
too. 

The stability of a crystalline phase is largely determined by the mini- 
mization of the Gibbs free energy [19]. We have compared the phase-wise free 
energy and the outcome is shown in Figure 9. We note from Fig. 9(a) that at 
300 K the free energy plot of the phase I joins smoothly at 6 GPa to that of 
phase Ha, which is consistent with the nature of second order phase transi- 
tion. Beyond 8 GPa, phase III has a lower free energy indicating the greater 
stability of this phase and also the onset of a first order phase transition from 
phase Ha to phase III. However, the transition pressure as obtained through 
the MD simulations [8] is higher than 8 GPa due to hysteresis. We observe the 
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greater stability of phase III at high pressures arises primarily due to its lower 
volume, while the vibrational energy remains lower and the entropy becomes 
higher in phase III providing an additional stability. From the plot of differ- 
ences in vibrational Gibbs free energy (AG) at 300 K as a function of pressure 
(Fig. 9(b)), we find that G for phase III is lower (above 8 GPa) than that 
for phase Ha. It is also observed (Fig. 9(c)) that AG between phase Ha and 
phase III structures decreases with temperature. Another point to note from 
the plots of relative free energy vs. temperature (at 8 GPa) is that in the very 
low temperature region (say, below 100 K), free energy difference shows an 
anomalous behaviour. It may be attributed to the fact that there is a greater 
density of low frequency modes in phase Ha than in phase III (due to the soft 
phonon modes in the former phase). Table 1 supports this view, where several 
vibrational modes have been compared between phase Ha and phase III. At 
low temperatures, these low frequencies are populated, thereby lowering the 
free energy. 

The above comparison of the Gibbs free energies in various phases pro- 
vides the phase diagram for the first order phase transitions. However, the 
second order transitions and also kinetic effects such as hysteresis in the first 
order transitions are better illustrated through molecular dynamics simula- 
tions. A qualitative impression of the phase transitions in LiYF4 is given in 
Fig. 10 with increasing pressure and temperature (see caption of Fig. 10). 
It stands our earlier observation [8] that phase II possesses two different space 
groups (viz. phase Ha: P1\jc ; phase lib: 12/ a) existing in two different ranges 
of temperaure. The low temperature phase Ha involves small displacements 
of some of the F atoms from their ideal positions of body-centered symmetry 
in phase lib. The transition Ha — > III occurs at a much higher pressure of 14 
GPa than the equilibrium phase boundary at 8 GPa due to hysteresis at the 
time scale of the simulation and is found irreversible on release of pressure. 



4 Conclusions 



We have demonstrated in this work that a lattice dynamics- molecular 
dynamics study can eventually lead one to calculate the vibrational free energy 
and other essential thermodynamic functions of a rare-earth ternary halide 
in various high-pressure phases. Effort has been made to shed light on the 
understanding of some important physical phenomena associated with the 
phase equilibria. Further, our results have the potential ingredients to help 
analyze the experimental inelastic scattering data, if high-pressure (along with 
low temperature) phonon measurements are to be carried out in future. 
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Table 1 

Comparison of calculated zone-center phonon frequencies between phase IIa(P= 13 
GPa) and phase III(P=18 GPa) of LiYF 4 at T=0 K. 

Long wavelength optical phonon modes (cm -1 ) 

Ag, Bg A u B u 



Phase Ha 


T ~\ 1 TTT 

Phase III 


Phase Ha 


I -\ l TTT 

Phase III 


Phase Ha 


T~»l TTT 

Phase III 


Phase Ha 


T~»1 TTT 

Phase III 


18 


92 


115 


127 














95 


107 


162 


151 


80 


82 








155 


151 


193 


165 


87 


128 


115 


142 


166 


175 


217 


200 


185 


156 


152 


183 


1 78 
1(5 


91 A 


997 


9^n 


iyo 


9nfi 

ZUD 


1 QQ 


1 

iyo 


215 


225 


260 


277 


266 


268 


220 


205 


286 


271 


293 


314 


293 


298 


277 


230 


303 


313 


316 


316 


310 


299 


304 


259 


324 


317 


336 


398 


343 


352 


310 


278 


334 


363 


369 


404 


358 


356 


315 


337 


371 


377 


384 


417 


408 


382 


345 


368 


413 


410 


438 


431 


434 


436 


370 


383 


429 


436 


461 


447 


455 


463 


388 


410 


435 


443 


496 


476 


463 


473 


437 


435 


473 


463 


507 


478 


495 


497 


477 


495 


506 


481 


531 


509 


514 


539 


494 


526 


546 


562 


577 


614 


578 


623 


553 


552 
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Table 2 

Calculated anisotropic thermal parameters (in units of 10~ 4 A 2 ) for the constituent 
atoms at various temperatures of 20 and 300 K associated with the three high- 
pressure phases of LiYF4. It may be noted that in the scheelite phase (74i/a, Z=4), 
all the F atoms are symmetrically equivalent. See text for the labeling of x, y, z 
directions. 







Phase I (P= 


=4 GPa) 


Phase Ha (P= 


T3 GPa) 


Phase III (P= 


-18 GPa) 


bp. 


Temp. 




















(K) 


U xx ^yy 


u 2Z 


^xx 


^yy 


u 2Z 




^yy 


v zz 


Li 


20 


63 63 


88 


54 


72 


60 


58 


64 


85 




300 


137 137 


224 


98 


144 


109 


95 


114 


189 


Y 


20 


17 17 


14 


13 


12 


17 


11 


12 


13 




300 


107 107 


71 


63 


52 


94 


39 


42 


51 


F(l) 


20 


41 41 


40 


36 


36 


36 


29 


28 


42 




300 


154 154 


124 


108 


103 


112 


70 


130 


106 


F(2) 


20 






36 


38 


42 


33 


30 


34 




300 






112 


113 


162 


88 


88 


106 


F(3) 


20 






36 


36 


38 


31 


36 


31 




300 






110 


108 


121 


78 


78 


106 


F(4) 


20 






35 


36 


35 


34 


45 


36 




300 






103 


100 


113 


93 


127 


106 
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Fig. 1. Crystal structures of LiYF4 belonging to (a) phase I, (b) phase Ila, (c) phase 
lib and (d) phase III. 

Fig. 2. Comparison of phonon dispersion relations along the same high symmetry 
direction of A in the three phases of LIYF4 at respectively 0, 13 and 18 GPa. Lines 
refer to our calculations while symbols represent the experimental data [2]. 

Fig. 3. (a) Calculated phonon density of states, g(E), in the three phases of LiYF4 
at respectively 4, 13 and 18 GPa; (b) comparison of g(E) in the two high-pressure 
phases at the same pressure of 13 GPa. 

Fig. 4. Calculated partial density of states for the constituent atoms with contri- 
butions from polarizations along three orthogonal directions in the three phases of 
L1YF4 at respectively 4, 13 and 18 GPa. Since phase I is of tetragonal symmetry, the 
contributions along x and y directions are identical for this phase. The relationship 
between the orthogonal directions and the crystallographic axes is given in section 
IIC. 

Fig. 5. Calculated (a) heat capacity (Cp) and (b) Debye temperature^^) as a func- 
tion of temperature for the three phases of LIYF4 at respectively 4, 13 and 18 GPa. 
The inset depicts the difference of Cp(T) and Cy(T) over the same temperature 
region. 

Fig. 6. Calculated average mode Griineisen parameter (7$) as a function of phonon 
energy for the three phases of LIYF4 at respectively 4, 13 and 18 GPa. 

Fig. 7. Calculated volume thermal expansion coefficient (ay) as a function of tem- 
perature for the three phases of LIYF4 at respectively 4, 13 and 18 GPa. (b), (c) 
and (d) demonstrate how phonons of different energies (in phase I, phase Ila and 
phase III respectively) contribute to ay at 20 and 300 K. 

Fig. 8. Calculated thermal mean squared displacements (on a semi-logarthmic scale) 
of the constituent atoms for (a) phase I, (b) phase Ila and (c) phase III due to 
phonons of various energies at 300 K. Insets in (a), (b) and (c) refer to the same but 
on an expanded scale up to 10 meV for the respective phases of LIYF4 at 4, 13 and 
18 GPa. It may be noted that in phases Ila and III, there are four symmetrically 
different F atoms, however, only the average contributions of all the F atoms are 
shown here. 

Fig. 9. (a) Comparison of the Gibbs free energy per atom as a function of pressure in 
all the three phases of LIYF4, as obtained from the lattice dynamical calculations 
at T = 300 K. The variation in Gibbs energy difference (AG) is plotted in (b) 
with pressure(P) and in (c) with temperature(T) giving a qualitative impression of 
the first order phase transition that takes place in switching from phase Ila to a 
dynamically more favourable phase III. Inset in (c) shows AG vs. T up to 200 K 
on an expanded scale. 

Fig. 10. Results of phase transitions as observed in MD simulations runs with in- 
creasing pressures at several constant temperatures (phase I — ► phase II(a/b) — > 
phase III) and also with increasing temperatures at several constant pressures (phase 
Ila — > phase lib). The transitions among the phase I, Ila and lib are found to be re- 
versible while phase III retains on decompression to P=0. The dashed line indicates 
the phase boundary between phase II and phase III as determined by quasiharmonic 
free energy calculations. 13 



